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We report measurements of the branching fractions for 0(3770) — * D°D°,D + D~, DD and res- 
onance parameters of ■0(3770) and tp(2S). By analyzing the line-shapes of the cross sections for 
inclusive hadron, D°D° and D + D~ event production in the range from 3.660 GeV to 3.872 GeV 
covering both 0(2/57) and 0(3770) resonances, we extract the branching fractions for 0(3770) de- 
cay into D°D° and D + D~ respectively to be 5(0(3770) -> D°D°) = (46.7 ± 4.7 ± 2.3)% and 
5(0(3770) D+D~) = (36.9 ± 3.7 ± 2.8)%, which give 5(^(3770) —*DD) = (83.6 ± 7.3 ± 4.2)% 
and non-ZXD branching fraction of 0(3770) to be 5(0(3770) -> non-DD) = (16.4±7.3±4.2)%. We 
meanwhile obtain the resonance parameters of 0(3770) and 0(2/57) to be M^( 3770 ) = 3772.2 ±0.7±0.3 
MeV, r^ 3770) = 26.9 ± 2.4 ± 0.3 MeV and F; e (3770) = 251 ± 26 ± 11 eV; M^ 2 s) = 3685.5 ± 0.0 ± 0.3 
MeV, r^ ( 2S) = 331±58±2 keV and Y^ (2S) = 2.330±0.036±0.110 keV; as well as the R value for the 
light hadron production directly through one photon annihilation to be R u ds = 2.262 ±0.054 ±0.109 
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PACS numbers: 



The 0(3770) resonance was discovered by MARK- 
I Collaboration in analysis of e + e~ annihilation into 
hadrons 29 years ago 0. Since its mass is above open 
charm-pair threshold and its width is two orders of mag- 
nitude larger than that of the ip(2S), it is thought to 
decay almost entirely to pure DD 2J. However, there 
were historically large discrepancies between the ob- 
served cross sections for DD and "0(3770) production 
in the e + e _ annihilation. The observed cross section 
°1a(3770) f° r 0(3770) production can be obtained based 
on 0(3770) resonance parameters and radiative cor- 
rections jj], which yield ^A 770 \ — 8.12 ± 1.56 nb at 
the center-of-mass (cm.) energy E c 



and 



_obs 

>(3770) 



7.53 ± 1.44 nb at E c 



= 3.7699 GeV 
= 3.773 GeV. 



There are three absolute measurements 00ll3 °f the 
observed cross sections for DD production in the e + e~ 
annihilation nearby the peak of 0(3770), which are sum- 
marized in table |U With the observed cross section 



T obs 
DD 



5.0 ± 0.6 nb for DD production at E c 



3.768 



GeV measured by MARK-III Collaboration, one histor- 
ically found that about 38% of 0(3770) does not de- 
cay into DD final states. Scaling the cr^g, measured 
at Epn = 3.768 GeV by using PDG 0(3770) parame- 
ters |3j yields the expected-observed cross section to be 
^scaled = 4j ±0 .6 nb at E cm = 3.773 GeV. The weighted 
average of the three cross sections listed in the second row 
of table Q] is = 6.25 ± 0.15 nb. Comparing the two 
observed cross sections for 0(3770) and DD production 
one finds that (17.0 ± 16.0)% of 0(3770) does not de- 
cay into DD, where the large error is dominated by the 
uncertainties in the resonance parameters of ■0(3770). 

To understand the discrepancy between the cr^g 770 ^ 
and o"^?, one has to more precisely measure both the 
0(3770) resonance parameters and the DD cross section, 
search for some 11011-DD decays of "0(3770) and directly 
measure the branching fraction for 0(3770) — * DD. To 
reduce some possible systematic shifts in comparison of 
the two cross sections from different experiments due to 
normalizations, one had better perform to simultaneously 
measure the "0(3770) resonance parameters and the DD 
cross section with a same data set. A better way to 
measure the inclusive branching fraction for 0(3770) — > 
DD is simultaneously to analyze the energy dependent 
cross sections for the inclusive hadron, D°D° and D + D~ 
event production in the energy range covering both 0(2,5) 
and 0(3770) resonances. In this way one can also more 
accurately measure the resonance parameters of both the 
0(3770) and 0(2$), since they are correlated each other 
in analysis of the cross section scan data. 

In this Letter, we present a line-shape analysis of the 
inclusive hadron, D°D° and D + D~ event production, by 
which we extract the branching fractions for 0(3770) — > 



TABLE I: The observed DD cross sections o"^,|, which were 
directly measured by MARK-III, CLEO and BES Collabora- 
tions, where scalcd means that the observed cross section is 
scaled from the measured value at 3.768 GeV to that at 3.773 
GeV based on PDG V0377O) resonance parameters. 

g cm (GeV) MARK-III (nb) CLEO (nb) BES (nb) 



3.768 
3.773 



5.0 ±0.6 
4.7 ± 0.6 scal< 



6.39 ±0.16 5.93 ±0.58 



for the first time, and meanwhile measure the resonance 
parameters of 0(3770) and ip(2S) with improved preci- 
sion on 0(3770) resonance parameters and with preci- 
sion in measuring leptonic width of 0(25*) comparable 
to the current PDG world averages. From this analysis 
we also measure the R value for light hadron production 
directly through one photon annihilation in the energy 
region from 3.660 GeV to 3.872 GeV. The analysis is 
based on the data taken with the BES-II detector [j| at 
the BEPC Collider during the time period from March 
to April, 2003. 

The observed hadronic cross section is determined by 
the relation 



obs 
^had 



/yobs 



r tri g : 

L Chad £>, 



(1) 



c had 



where is the number of the observed hadronic 

events, L is the integrated luminosity, Chad is the effi- 
ciency for detection of the inclusive hadronic events and 
e had i s the trigger efficiency for collecting the hadronic 
events in on-line data acquisition. 

The hadronic events are required to have more than 
2 good charged tracks, each of which is required to 
satisfy the following selection criteria: (1) the charged 
track must be with a good helix fit and the number of 
dE/dx hits per charged track is required to be greater 
than 14; (2) the point of the closest approach to the 
beam line must have radius < 2 cm; (3) |cos6*| < 0.84, 
where 9 is the polar angle of the charged track; (4) 
2.0 ns < Ttof < T p + 2.0 ns, where Ttof is the 
time-of-fiight of the charged particle, and T p is the ex- 
pected timc-of-flight of proton with a given momentum; 

(5) the charged track must not be identified as a muon; 

(6) p < E b + 0.1 x E b x ^/(l + El), where p is the charged 
track momentum and E\, is the beam energy in GeV; (7) 
for the charged track, the energy deposited in the BSC 
should be less than 1.0 GeV. In addition, the total en- 
ergy deposited by an event in the BSC should be greater 
than 28% of the beam energy. Furthermore, the selected 
tracks must not all point into the same hemisphere in the 
z direction. 

Some beam-gas associated background events can also 
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FIG. 1: (a) the distribution of the averaged z of the events 
satisfying the hadronic event selection criteria, where the his- 
togram shows the events from the data, and the curves give 
the best fit to the z distribution; (b) the efficiencies for detec- 
tion of the inclusive hadronic events vs the nominal center-of- 
mass energies. 



ground events are produced at random z positions in 
the beam pipe, while the hadronic events are produced 
around z — 0, this characteristic can be used to distin- 
guish the hadronic events from the beam-gas associated 
background events. Fig. Qfa) shows the distribution of 
the averaged z of the accepted events which satisfy the se- 
lection criteria. Using a Gaussian function plus a second 
order polynomial to fit the averaged z distribution of the 
events, we obtain the number, iVJ^J, of the candidates for 
hadronic events. This number of candidates for hadronic 
events contains some contaminations from some back- 
ground events such as e + e~ — > t + t~ , e + e~ — > ("f)e + e~ , 
e + e~ — > (7)/i + /i _ , and two-photon processes. The num- 
ber of the background events, Nb, due to these processes 
can be estimated by means of Monte Carlo simulation. 
Subtracting from N?^ yields the number of the ob- 
served hadronic events, iVJ~|. The systematic uncer- 
tainty in measuring the produced hadronic events due 
to the hadronic event selection criteria is estimated to be 
about - 2.5%. 

The integrated luminosities of the data sets are deter- 
mined using large-angle Bhabha scattering events, which 
satisfy the following selection criteria: (1) two charged 
tracks with total charge zero, for each track, the point 



radius < 1.5 cm and \z\ < 15 cm where \z\ is measured 
along the beam line from the nominal beam crossing 
point; (2) each track is required to satisfy | cos#| < 0.7, 
where 9 is the polar angle of the charged track, (3) it is re- 
quired that the energy deposited for each charged track 
in the BSC be greater than 1.1 GeV and at least the 
magnitude of one charged track momentum be greater 
than 0.9E^, (4) no track goes through the rib regions of 
the BSC. The systematic uncertainty in the measured lu- 
minosities arises mainly from the difference between the 
data and Monte Carlo simulation. The total uncertainty 
is estimated to be ~ 1.8%. 

The detection efficiency for hadronic events is deter- 
mined via a special Monte Carlo generator m which 
the radiative corrections to a 2 order are taken into ac- 
count. These generated events are simulated with the 
GEANT3-based Monte Carlo simulation package [llj| . 
The systematic uncertainty in the efficiencies due to the 
generator is estimated to be ~ 2.0% (~ 0.7%) for re- 
construction of the hadronic events from -0(3770) and 
ip(2S) decays (from continuum hadrons). Fig.QJb) shows 
the Monte Carlo efficiencies for detection of the hadronic 
events produced at the different nominal center-of-mass 
energies. These efficiencies are used in the measurements 
of the observed inclusive hadronic cross sections at each 
of the energy points. 

The trigger efficiencies are obtained by comparing the 
responses to different trigger requirements in the data 
taken at 3.097 GeV during the cross section scan ex- 
periment. The efficiencies are measured to be etrig = 
(100.0t^)% for both the e+e" -> (j)e + e~ and e+e" -> 
hadrons events. 

The observed cross section for D°D° (or D + D~) pro- 
duction is determined based on the number N n ° (or 
N n + ) of the reconstructed D° (or D + ) events by 

tag 

N D o (or N n + ) 

0bs_ _ -"tag ^ -L'tag ' fr> s 

° DODO (or D+D-) ~ 2 xLxBxe ' 1 ' 

where L is the integrated luminosity of the data set used 
in the analysis, B is the branching fraction for the de- 
cay mode in question, and e is the efficiency determined 
from Monte Carlo simulation for reconstruction of this 
decay mode. The observed numbers, N D o (or N D + ), 

of the singly tagged D° (or D + ) are obtained by analyz- 
ing the invariant mass spectra of K^n^ and ir^ ir^ 
(or K^ 1 ^^) as discussed in detail in Ref. [T2|. As an 
example, Fig. shows the distributions of the invariant 
masses of K~ir + combinations observed from 14 data sets 
collected at different energy points in ^(3770) resonance 
region. 

Using the methods discussed above and in Ref. (l2^ . 
we obtain the observed cross sections for both the inclu- 
sive hadron and DD event production at the energies at 
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FIG. 2: The distributions of the invariant masses of the K~tt 
combinations from the data sets taken at different nominal 
cm. energy. 



cross sections (points with error) for inclusive hadronic 
event production, while Fig. EJb) and 0Jc) respectively 
display the observed cross sections (points with error) 
for D°D° and D + D~ production, where the error bars 
represent the combined statistical and point-to-point sys- 
tematic uncertainties including the statistical uncertain- 
ties in the luminosity and the efficiencies for detection of 
the hadronic events and Bhabha events. 

The determination of the branching fractions for 
V>(3770) -> D°D° 1 D+D~, and DD is accomplished 
by simultaneously fitting the observed cross sections for 
ip(2S), ^(3770), D°D° and D+D~ to functions that de- 
scribe the combined ip(2S), i/'(3770) resonance shapes, 
the tail of J/ip resonance and the non-resonant hadronic 
background, as well as the partial ■0(3770) resonance 
shapes for 0(3770) -> D°D° and 0(3770) -> D+D~. 
The functions are corrected for the radiative correc- 
tions 30. 
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FIG. 3: The hadronic cross sections versus the nominal center- 
of-mass energies, where (a) the points with error show the 
observed cross sections, the solid line shows the fit to the cross 
sections and the dashed line represents the contributions from 
J/ip, ^(2S) and continuum hadron production; where (b) the 
points with error show the observed cross sections and the 
line 1 gives the fit to the data in ij>(2S) resonance region, 
while the squares with error show the observed cross sections 
for the continuum hadron production, which are corrected for 
the radiative corrections; the line 2 gives the fit to the cross 
sections (see text). 



Breit-Wigner function 

B , l2nT ee T h 

a (s ' ~ (7 - m 2 ) 2 + (r tot M) 2 ' (3) 

to describe their production, where s is the squared ac- 
tual energy which produces the hadronic events, M and 
ptot are respectively the masses and total widths of the 
resonances, and T ee and T h are the partial widths to 
e + e~ channel and to the inclusive hadronic final states, 
respectively. Assuming that there are no other new struc- 
tures and effects except the 0(3770) resonance and the 
continuum hadron production in the energy region from 
3.70 GeV to 3.86 GeV, we use the pure p-wave Born 
order Breit-Wigner function with an energy-dependent 
total width to describe the ^(3770) production and the 
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FIG. 4: The observed cross sections versus the nomi- 
nal center-of-mass energies, where (a) shows the inclusive 
hadronic event production, (b) and (c) show the D°D° and 
D + D~ event production, respectively; the points with error 
are the data, while the lines are the fits to the data. 



decays. The -0(3770) resonance shape is taken as 



)(*) 



1 9-!rV ee F tot (1 } 
iZ7ri i/>(3770r <£(3770)V 6 > 



^ 3my " ' ~ ^(3770)) 2 + [^(3770)r^ 770) ( S ')] 2 ' 

( 4 ) 

while the D°D°, D + D (or DD) resonances shapes are 
taken as 



127rr T(3770) r -D-D( S ) 



(«' - M |(3770)) 2 + [^(3770 ) r t ; ( t 3770) ( S ')]' 



(5) 

where M^( 3770 ) and r^ e (3 77 o) are ^ e mass an d leptonic 
width of the 0(3770) resonance, respectively; is the 
partial width of 0(3770) decay into DD; Fffi 3770 Js ) and 

Tdd( s ) are chosen to be energy dependent, which are 
defined as 



r$(3770)( s ) — r D o D o(s ) + T D+D - (s ) + r noa __ D _ D (s 



where 



'd»d»( s ) — r 0d°d° , in n 



(6) 



{p DO fl + {rp%)\ 

Boo, (7) 



(p J+ ) 3 l + (r P ° D+ ) 2 



where pj, and p£> are respectively the momenta of the 
D mesons produced at the peak of 0(3770) and at the 
actual cm. energy W ; To is the total width of the 
0(3770) at its peak, B Q0 = £(0(3770) -> D°D Q ) and 
B+- = 5(0(3770) -> D+D~) are the branching frac- 
tions for 0(3770) -> Z?°Z) and V(3770) -► £>+£>-, re- 
spectively, which are the fitted parameters, and r is the 
interaction radius of the cc, which is left free in the fit; 
@d°d° an d ®d+ D~ are the step functions to account for 
the thresholds of the D°D° and D + D~ production, re- 
spectively 

The non-resonant hadronic background shape is taken 

as 



nrsnt 
Oh 



ds"G(s, s") r ^y^^ FM 



/() 

+f D D 



|i-n( S ') 



r PD° \3 a , (Pd+\3q 



D 



(10) 



with s' = s(l — x), where a; is a parameter related to the 
total energy of the emitted photons and y/s is the nom- 
inal cm. energy, F(x,s) is the sampling function U 
1/1 1 — II(s(l — x))\ is the vacuum polarization correc- 
tion function []Jfl including the contributions from all 
1 resonances, the QED continuum hadron spectrum 
as well as the contributions from the lepton pairs (e + e~, 
and t + t~) |10J; crf[ + (s 1 ) is the Born cross sec- 
tion for e + e~ — > /i + /J~, -Erjo and £x>+ are respectively 
the energies of D° and Z? + mesons produced at the ac- 
tual energy yfs 1 ', /p^ is a parameter to be fitted, and 
Ruds(s') is the R value for the light hadron production 
through one photon annihilation directly. In the fit we 
take the R u ds(s') as a constant in the energy region and 
left it free; we fix the J/0 resonance parameters at the 
values given by PDG 0- We also consider the effects of 
the BEPC energy spread on the observed cross sections 
in the fit. G(s, s") in Eq. (10) is the Gaussan function 
to describe the cm. energy distribution of the BEPC 
machine. 

The curves in Fig.[3]and Fig.0]show the fits to the data. 
Fig. EJa) shows the observed cross sections with the fit to 
the data, where the points with error show the observed 
cross sections and the error is combined from statistical 
and point-to-point systematic uncertainties arising from 
the statistical uncertainties in the efficiencies for detec- 
tion of the hadronic events and Bhabha events; the solid 
line shows the fit to the cross sections and the dashed 
line represents the contributions from J/0, ip("2S) and 
continuum hadron production. To examine the contri- 
bution from the vacuum polarization corrections to the 
Born hadronic cross section due to one photon annihila- 
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-0(3770) as well as J/-0 from the observed cross sections 
to yield the expected cross sections of the continuum 
hadron production corrected with the radiative effects, 
which is given by Eq. (10). The squares with error in 
Fig-EIb) show the yielded-expected cross sections, where 
the errors are the originally absolute errors of the totally 
observed cross sections as shown in Fig.^a). The line 2 
in Fig. Gfb) shows the fit to the expected cross sections 
of the continuum hadron production corrected for the 
radiative effects. Fig. 0Ja) shows the observed cross sec- 
tions for the inclusive hadronic event production, where 
the contributions from J and ip{2S) radiative tails as 
well as the continuum hadron production are removed. 
Fig-Gib) and Fig.0fc) display the observed cross sections 
for D°D° and D + D~ production, respectively, where the 
points with error show the observed cross sections, while 
the lines give the fits to the data. The error is combined 
from the statistical and point-to-point systematic aris- 
ing from the statistical uncertainties in the efficiencies 
for detection of the singly tagged D events and Bhabha 
events. The fit gives x 2 / n - oi = 65.4/64 = 1.02. In the 
data reduction, the number of the singly tagged D° (or 
D + ) events are removed from the inclusive hadronic event 
samples before calculating the hadronic cross section and 
its error based on Eq.(l). These make the hadronic event 
samples and the singly tagged D samples be independent. 
The results from this fit are 



and 



5(0(3770) -> D°D°) = (46.7 ± 4.7 ± 2.3)% 



5(0(3770) -> D+D-) = (36.9 ±3.7 ±2. 



where the first error is statistical and second system- 
atic arising from uncanceled systematic uncertainties in 
the measured Chad (~ 2.8%) in (Tjjoqo (~ 4.1%), and in 
a D+D- (~ 7.0%). Considering the correlation between 
the two branching fractions obtained from the fit, we ob- 
tain the branching fraction for V>(3770) — > DD to be 

5(0(3770) -> DD) = (83.6 ± 7.3 ± 4.2)%, 

which results in the non-5Z) branching fraction to be 



5(0(3770) 



non 



DD) = (16.4 ±7.3 ±4.2)%. 



The fit gives the BEPC machine energy spread oe bepc = 
(1.343 ± 0.029) MeV. From the fit we obtain the R value 
for the light hadron production due to one photon an- 
nihilation in the energy region between 3.660 GeV and 
3.872 GeV to be 

R uds = 2.262 ± 0.054 ±0.109, 

where the first error is statistical and the second system- 
atic arising from the uncertainty in the measured cross 



(~ 0.2%). The fit also gives the resonance parameters of 
0(3770) and 0(25) to be M^rm) = 3772.2 ± 0.7 ± 0.3 
Me V, T% m = To = 26.9 ± 2.4 ± 0.3 MeV, r£ (3770) = 
251 ± 26 ± 11 eV; M^s) = 3685.5 ± 0.0 ± 0.3 MeV, 
r ^°(2S) = 331±58±2 keV, r^ e (2S) = 2.330±0.036±0.110 
keV, where the first error is statistical and the second sys- 
tematic arising from the uncertainties in the measured 
Chad (~ 4.4%), in a D of>o (~ 4.5%), and in a D+D - (~ 
7.4%). The measured mass difference between 0(3770) 
and 0(25) is AM = 86.6±0.7 MeV. The measured values 
of the -0(3770) resonance parameters yield the cross sec- 
tion for 0(3770) production at its peak to be cr^ r ( 3770 ) — 
(9.63 ± 0.66 ± 0.35) nb, corresponding the observed cross 
(6.94 ± 0.48 ± 0.28) nb, which is con- 
;8.12 ± 1.56) nb at 
0(3770) peak obtained based on the PDG 0(3770) reso- 
nance parameters. 

The measured branching fractions yield the ratio of 
the partial widths T DO£) o/T D+D - = 1.27 ± 0.12 ± 0.08, 
which agrees with the prediction of 1.36 by Eichten et 
al. ^ij and is in good agreement with 1.41 ± 0.23 ± 0.11 
measured by BES 15]. The measured widths of the 
resonances yield the leptonic branching fractions to be 
55(0(3770) -> e+e-) = (0.93 ± 0.06 ± 0.03) x 10~ 5 and 
BF(ip(2S) -► e+e") = (0.704 ± 0.122 ± 0.033)%. 

We find that the continuum background shape affect 
the measured total and leptonic widths of the resonances 
from the line-shape analysis. If we take the non-resonant 
hadronic background shape as 



section <J° h { l 770) 

sistent within error with c^g 770 ^ 



nrsnt 



h a.. 



+ 



V) 

fan [iw^f 1 



DODO + 



B 



-(»'), 



(11) 



in fitting the data (where h is a parameter in the fit), 

290 ± 59 ± 5 keV, Tf, 



we would obtain r^g) 



- V(2S) 



2.378 ± 0.036 ± 0.103 keV, T^ 377Q) = 27.3 ± 2.5 ± 1.1 
MeV and r^f (3770) = 256 ± 27 ± 13 eV, and almost un- 
changed the measured masses of the resonances. This fit 
yields x 2 / n -o.f = 75.3/64 = 1.18. These indicate that 
the vacuum polarization corrections to the Born order 
cross sections for the continuum hadron production can 
not be ignored in precisely measuring the resonance pa- 
rameters of the narrow resonances like J/tp and ip(2S) 
as well as T(15) etc. in e + e~ cross section scan exper- 
iments. Ignoring the effects of the vacuum polarization 
corrections on the continuum hadron production cross 
sections in analysis of the cross section scan data taken 
in the tp(2S) resonance region would decrease the ip(2S) 
total width by about 40 keV. 

In summary, we measured the branching fractions for 
0(3770) -► D°D°,D+D- and V(3770) -> non - DD for 
the first time and measured the resonance parameters of 
0(3770) and 0(2S*) with improved precision on -0(3770) 
resonance parameter and with a comparable precision to 
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same data samples, we also measured the R value for 
the continuum light hadron production with a precision 
of about ±5%. From this analysis we directly observed 
the effects of the vacuum polarization on the observed 
cross sections of the continuum hadron production in the 
neighborhood of the ip(2S) resonance. 
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